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Abstract. Open charm spectra at mid-rapidity from direct reconstruction of 
D° — > Kir and indirect electron/positron measurements via charm semileptonic 
decay in p+p and d+Au collisions at ^/s NN = 200 GeV are reported. The 
combined spectra cover open charm transverse momentum range 0.1^j»t^6 GeV/c. 
The electron spectra show approximate binary scaling between p+p and d+Au 
collisions. From these two independent analyses, the differential cross section per 
nucleon-nucleon interaction at mid-rapidity for open charm production at RHIC 
is dcr^ Ar /d?/=0.31±0.04(stat.)±0.08(syst.) mb. The next-to-leading-order pQCD 
calculation underpredicts the charm yields. 

The production and spectra of hadrons with heavy flavor are sensitive to initial 
conditions and the later stage dynamical evolution in high energy nuclear collisions, and 
may be less affected by the non-perturbative complication in theoretical calculations pQ. 
Charm production has been proposed as a sensitive measurement of parton distribution 
function in nucleon and the nuclear shadowing effect by systematically studying p+p, 
and p+A collisions 2\. The relatively reduced energy loss of heavy quark traversing 
a Quark-Gluon Plasma will help us distinguish the medium in which the jet loses its 
energy A possible enhancement of charmonia (J/^f) production can be present at 
RHIC energies jl] due to the coalescence of the copiously produced charm quarks. 

Direct and indirect charm production has been measured in hadron-hadron 
collisions El El • However, no experimental results are available at v /s NN =200GeV. 
Theoretical predictions for this energy region differ significantly 8J. Therefore precise 
experimental measurement of the baseline charm cross sections at RHIC is necessary. In 
this paper, we report the recent STAR results on the absolute open charm cross section 
measurements from direct charmed hadron D° reconstruction in d+Au collisions and 
electrons from charm semileptonic decay in both p+p and d+Au collisions at 200 GeV. 

1. Analysis methods and results 

At STAR, the hadronic channel of D° — > K~tt + (+c.c.) with branching ratio of 3.83% 
was measure by the Time Projection Chamber (TPC) [HI EH] in d+Au collisions. The 
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Figure 1. (a) Invariant mass distributions of oppositely charged Kit pairs from d+Au 
collisions, (b) dE/dx in TPC versusp with the TOFr velocity cut \l/(3-l\ < 0.03. The 
insert shows dE/dx distribution for 1 < p < 1.5 GeV/c. (c) The e + e~ pair invariant 
mass square distribution from different sources in PYTHIA |15| simulation. 



D° signal was reconstructed using the event-mixing technique which is the same as 
the one used to identify and measure the K°* resonance [llj . After the mixed-event 
background subtraction, the D° signal in p? < 3 GeV/c and \y\ < 1 from 15.7 million 
d+Au minimum bias events is shown as filled symbols in Fig. 1(a). A linear function was 
used to reproduce the residual background. After the two-step background subtraction, 
a gaussian function was used to obtain the D° signal shown as open symbols. 

In the 2002-2003 run, a prototype multi-gap resistive plate chamber time-of-flight 
system (TOFr) was installed with the coverage — 1 < r\ < in pseudorapidity. In 
addition to its capability of hadron identification ^2j, electrons could be identified at 
low momentum (px < 3 GeV/c) by the combination of velocity {(3) from TOFr |Hj 
and the particle ionization energy loss (dE/dx) from TPC 0. Fig. 1(b) shows that 
the electrons are clearly identified as a separate band in the dE/dx versus momentum 
(p) with a selection on /3 at |l//3 — 1| < 0.03. At higher p T (2-4 GeV/c), negative 
electrons were also identified directly by TPC since hadrons have lower dE/dx due to 
the relativisitic rise of electron dE/dx. Detector acceptance and efficiency corrections 
were determined from the embedding data [T^]. Total inclusive electron spectra from 
p+p and d+Au collisions are shown as symbols in Fig. 2. 

7 conversions and tt° Dalitz decays are the dominant photonic sources for electron 
background. Such background was mainly constrained at small pair invariant mass and 
small opening angle. To measure the background, pair invariant mass and opening angle 
distributions were constructed by first selecting an electron from TOFr and then looping 
through oppositely charged tracks reconstructed in the TPC Fig. 1(c) shows an 
example of the pair invariant mass. Simulations using both HIJING ^3] and PYTHIA 
[TH] . with GEANT to describe the detector, showed ~60% efficiency of such background 
reconstruction for p T > lGeV/c, with negligible p T dependence. This value was used 
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Figure 2. Electron distributions from p+p (left) and d+Au collisions (right), 
respectively. The symbols are electrons identified via TOFr combined with TPC. 



to correct for the photonic electron spectra in data. In this method, more than 95% 
background has been measured, shown as thick solid lines in Fig. 2. Other contributions 
(< 5%) from decays of r],uj,p,<j) and K were determined from simulations. The overall 
uncertainty of the background is ~20%. 

The non-photonic electron spectra, obtained by subtracting the previously 
described photonic background from the inclusive spectra, are shown as symbols in 
Fig. 3 (left) for both p+p (triangles) and d+Au (circles) collisions. The D° invariant 
yields d 2 N / 2npTdpTdy are calculated and also shown in Fig. 3 (left) as squared symbols 
as a function of pt- The mid-rapidity yield dN/dy was extracted with an exponential fit 
to the invariant yield as a function of transverse mass m?, which is 0.028 ±0.004 ±0.007. 
Several model studies PJ showed that semi-leptonic decays from open charm are the 
dominant non-photonic electrons at 1 < pj- < 4 GeV/c. We performed a fit with the 
combined results of D° and electron distributions in d+Au collisions, assuming that the 
D° spectrum follows the power law up to 6 GeV/c. From the fit, we got the yield of 
D° at mid-rapidity dN/dy = 0.030 ± 0.004 ± 0.008 and dcrg N /dy= 0.31 ± 0.04 ± 0.08 
mb. Thus the total charm-pair cross section a^ 5 N = 1.44 ± 0.20 ± 0.44 mb, where a 
factor of 4.7 ± 0.7, estimated by the model simulations IBHH3, was used to convert the 
da/dy at mid-rapidity to total cross section. The systematic error is dominated by the 
background subtraction, p? coverage for each measurement and overall normalization. 
The nuclear modification factor was obtained by taking the ratio of the electron spectra 
in Fig. 3 (left) scaled with the underlying nucleon-nucleon binary collisions. It was 
measured to be 1.2 ± 0.2 ± 0.3 and is consistent with binary scaling within the errors. 

The beam energy dependence of the cross section is shown in Fig. 3 (right). 
The solid (MSEL=1) and dashed (MSEL=4) lines are PYTHIA calcualtions with and 
without high order processes such as flavor excitation etc., respectively. The dot-dashed 
line is the NLO calculation from [H], which underpredicts our result. Dotted line is the 
power-law fit of a^ 5 N oc (^fs) n to the data points. For the total charm production, the 
power n ~ 2.2±0.5, while n ~ 0.5(0.3) has been observed for pion (charged multiplicity) 
productions at lower energy 01 ^] . 
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Figure 3. (left) Reconstructed D° px distributions from d+Au collisions. Non- 
photonic electron pt distributions from p+p and d+Au collisions (open and filled 
symbols are electrons identified via TOFr combined with TPC and TPC only). Solid 
and dashed lines are the combined fit results for the D° and electron spectra in d+Au 
collisions. Dotted line is scaled down by the number of binary collisions (iV&i n =7.5) in 
d+Au collisions, (right) Total charm-pair cross sections versus the collision energy (in 
• v /s NN ). The result of our measurement is shown as star. 



2. Summary 

The charm cross section, d(T^ JV /dy=0.31±0.04(stat.)±0.08(syst.) mb from 200 GeV 
d+Au collisions has been measured at RHIC. The independent measurements of the 
reconstructed D° and single electrons from charm semileptonic decay are consistent. 
The NLO calculation underpredicts the total cross section at this energy. 
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